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' Abstract 

In this paper, we derive optimal upper and lower bounds on the dimension of the at- 
tractor Aw for scalar reaction-diffusion equations with a Wentzell (dynamic) boundary 
, condition. We are also interested in obtaining explicit bounds on the constants in- 

I volved in our asymptotic estimates, and to compare these bounds to previously known 

[-^ ' estimates for the dimension of the global attractor Ak^ K € {D,N,P}, of reaction- 

diffusion equations subject to Dirichlet, Neumann and periodic boundary conditions. 
I The explicit estimates we obtain show that the dimension of the global attractor Aw 

■ is of different order than the dimension of Ak, for each K £ {D,N,P} , in all space 

dimensions that are greater or equal than three. 



X 

^ ; 1 Introduction 

It is well-known that the long-time behaviour of solutions of partial differential equations 
arising in mathematical physics can, in many cases, be described in terms of global attractors 
of the associated semigroups (see |3l El |32l |16] and references therein). For a large class of 
equations of mathematical physics, including parabolic partial differential equations mod- 
elling reaction, diffusion and drift, hyperbolic type equations, and so on, the corresponding 
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attractor has finite Hausdorff and fractal dimensions. Tlius, tlie dynamics on tlie attractor 
liappens to be finite-dimensional, even though the system is governed by a set of partial dif- 
ferential equations. As the dimension of the attractor is indicative of the number of degrees 
of freedom needed to simulate a given dynamical system, it is then crucial to obtain more 
realistic estimates for its dimension in terms of observable physical quantities. 

Aside from some applied motivation, much of the mathematical interest nowadays is 
centered on the dynamics of boundary value problems with static boundary conditions of 
Dirichlet and Neumann-Robin type, or even periodic boundary conditions. The influence of 
these dissipative boundary conditions on a given model has only been recently investigated 
in connection with a class of react ion- diffusion systems. In [33], a first contribution is made 
to the understanding of this problem with a Robin boundary condition. In particular, it 
is shown, for a fixed nonlinearity, how the flow defined by the reaction-diffusion system 
depends on the interaction between diffusion u and another parameter 6 involved in the 
boundary condition (cf. also |34])- A classification of points in (z/, 6')-space, as structurally 
stable, or bifurcation points, for a one-dimensional scalar reaction-diffusion equation with a 
cubic nonlinearity is discussed in detail in [33]. Other studies on the influence of boundary 
conditions upon the solution structures of partial differential equations have also been done 
by other scientists. These studies have analyzed the detailed effect of boundary conditions on 
the structure of global attractors (see, e.g., [HI EQl |36l E] ) . If the equilibrium is nonhyperbolic 
and a bifurcation occurs, the structure of attractors may vary with respect to boundary 
conditions. This has been observed in the analysis of pattern formation in a ID reaction- 
diffusion system [8], in lattice systems [53], in the study of steady state bifurcations [30] 136]. 
and finally in [H], on mode-jumping of the von Karman equations. 

Although the global attractors of these systems will depend, for a given nonlinearity, on 
the choice of the boundary conditions, their finite dimension does generally not. This result 
can be easily formulated for a scalar reaction-diffusion equation, as follows. Consider the 
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parabolic partial differential equation 



dtu = uAu - f (u) + Xu + g, {x,t) e Q x {0, +00) , (1.1) 

where u = u{x,t) G M, 1^ C M", n > 1, is a bounded domain with sufficiently smooth 
boundary T, g = g (x), and u, A are positive constants. The function / : M — )■ M is assumed 
to be C^'^, that is, continuous and with a Lipschitz continuous first derivative, which satisfies, 
among other natural growth conditions (see, e.g., [6l Chapter II]), 

/ (y) > — c/, for all y G M, for some C/ > 0. 

We may ask that u satisfy either a Dirichlet {K = D) boundary condition or a Neumann 
{K = N) boundary condition, and even a periodic {K = P) boundary condition. It is well- 
known that equation (11. ip . supplemented with an appropriate initial condition, generates a 
semigroup {St} acting on a suitable phase space H. This semigroup possesses the global 
attractor Ak, which may depend on the choice of the boundary conditions, and Ak has 
finite fractal dimension for each K G {D,N,P}. In particular, the Haussdorf and fractal 
dimensions of Ak, for any K G {D, N, P}, satisfy the following upper and lower bounds: 

-j \n\ <dimHAK <dimFAK <ci(l + ^^^^\n\^^''] , (1.2) 

for some positive constants Cq, Ci that depend only on n, f and the shape of Q (see, e.g., [3l 
Chapter III]; cf. also [6], [Ml Chapter VI]). Here, \Q\ stands for the Lebesgue measure of 
Q. For a fixed domain Q, we observe that these estimates are sharp with respect to z/ — )■ 0^ 
(for each fixed A > 0), or sufficiently large A (for each fixed u > 0). Hence, these bounds 
for the dimension of Ak are of the same order for each K G {D, N, P}. These remarkable 
estimates also depend linearly on the "volume" of the spatial domain Q, which is consistent 
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with physical intuition. This property of the dimension of the attractor has not been proved 
for all equations, such as, the Kuramoto-Sivashinski equation. 

Our main goal in this paper is to investigate the dependance of the dimension of the global 
attractor for equation (11.11) subject to a completely new class of boundary conditions, which 
are sometimes dubbed as Wentzell boundary conditions, and which have some applications 
in probability theory, specifically, Markov processes. But what are they really? To put 
them into a context, let L be an elliptic differential operator of the second order (e.g., 
L = uA) with coefficients that are well-defined over Q. It is known that there is a one-to-one 
correspondence between (Co)-semigroups and Markov processes in Q which are homogeneous 
in time and satisfy the condition of Feller |9] (that is, the range of the resolvent operator 
coincides with a prescribed set). Thus, to each such Markov process there is a corresponding 
semigroup of operators 

TtV (x) = j v{y)P{t,x,dy), 
n 

where the Markov transition function P{t,x,B) satisfies P{t,x,B) > 0, for t > 0, x E fl 
and any Borel set B C Q. As a function of B, P {t, x, ■) is a probability measure. What are 
the most general boundary conditions which restrict the given operator L (more correctly, its 
closure) to the infinitesimal operator of a semigroup of positive contraction operators acting 
on C (^])? Wentzell ^ gave a partial answer to this question in higher space dimensions by 
finding a sufficiently large class of boundary conditions which involve differential operators 
on the boundary that are of the same order as the operator acting in Q. He discovered the 
following form of boundary conditions: 

Lu + vhd^u = 0, on r X (0, +oo) , (1.3) 

where n denotes the outward normal at F, 6 is a positive constant and d^u is the outward 
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co-normal derivative of u with respect to L. We refer also to the pioneering work of |10j . 
for generation theorems for L with Wentzell boundary conditions in one space dimension. 
Until the work of [H], the study of the operator L with Wentzell boundary conditions was 
usually confined to generation properties of this operator in the space C (Q) . In 2002, the 
authors in [11] have found a way to introduce the Wentzell boundary condition (11.31) in an 
L^-context, which led to the discovery of the natural space for these type of problems (see 
Section 2). The reader is referred to |11[27] for an extensive survey of these results and some 
history. 

For the homogeneous linear heat equation (II. ip (that is, f = g = X = 0), the Wentzell 
boundary condition (II. 3p is equivalent to a purely differential equation of the form 

dtu + ubdnu = 0, onrx(0,oo). (1.4) 

Thus, the main attraction here is that there is a dynamic element introduced into the bound- 
ary condition. The heat equation, supplemented by either (11.31) or (II. 4p . corresponds to the 
situation where there is a heat source (if 6 > 0) or sink (if 6 < 0) acting on the bound- 
ary r. Mathematically speaking, this kind of conditions (ll.4p arises due to the presence of 
additional boundary terms in the free energy, which must also account for the action of a 
source on F (see [22]). We refer the reader to p6^j (cf. also [16]), for an extensive derivation 
and physical interpretation of (II. 4p for (II. ip . For the nonlinear parabolic equation (II. ip . the 
boundary condition 01. 3p can be formally be transformed into a condition of the form 

dtu + vbdnU + f (u) — Xu = g, on F x (0, oo) . (1.5) 

Generally, one may replace / — A in (II. 5p by another arbitrary function h, satisfying suitable 
assumptions. With more sophisticated arguments, using techniques from semigroup theory, 
and a variation of parameter formula, it is possible to prove that the regularity of the solution 
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for (ll.ip .( !L5|) increases as /, Q and g become more regular (see Section 2). In particular, for 
(yf = 0, if f2 is a bounded C°° domain and / is a C°° function, regularity theory implies that 
the solution u (t) to f ll.ip .( lL5l) belongs (Q) , for all /c > and all positive times. At least 
in this case, the boundary condition {\1.3\i . for equation (11 .11) . is equivalent to the boundary 
condition (ll.Sp . However, in general, this may not be so if the solution, for the semilinear 
problem (II. ip and the Wentzell condition (II. 5p . is not smooth enough. Since we wish to 
treat the most general case, by imposing the least regularity assumptions on /, g and Q, we 
will devote our attention only to the study of (II. ip . subject to linear boundary conditions of 
the form (II. 4p . Our results below can be immediately extended to other classes of nonlinear 
Wentzell boundary conditions (see, e.g., [22] and references therein). Boundary conditions of 
the form (II. 5p arise for many known equations of mathematical physics. They are motivated 
by heat control problems formulated in the book of Duvaut and Lions [7] , problems in phase- 
transition phenomena [HI [151 113 EOl EH EU [251 E3 [M] (and their references), special flows in 
hydrodynamics [12l [22l [42l [39], Stefan problems [H [35l [41], models in climatology [40], and 
many others. The reader is referred to for a more complete list of references involving 
the application of such boundary conditions to real-world phenomena. 

By keeping our treatment of the boundary condition simple, we wish to prove that the 
problem (II. ID . (II. 4p generates a dynamical system on a suitable phase-space, possessing a 
finite dimensional global attractor Aw- Then, we establish that the Haussdorf and fractal 
dimensions of Aw satisfy the following upper and lower bounds: 

A / _|_ jSy \n-l 

Cw{n,T)u ) - ^^"^^ - ^^"^^ - + CwU r) J ' ^'-'^ 

for n > 2, and 

/ ^ \ 1/2 / c + A \ ^''^ 

^3(7^^?^) < dimn Aw < dimp Aw<cjl+ , „ , I , (1.7) 
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in one space dimension. The positive constants Cj, i = 1, •••,4, depend only on n, / and the 
shape of Q, while explicit estimates and formulas for Cw T) and Cd (Q) , respectively, 
are provided in the Appendix. We note again that, for a fixed domain Q, these estimates 
are sharp with respect to — )■ 0^ (for each fixed A > 0), and for sufficiently large X {if u > 
is fixed). We remark that the bounds we obtain in fll.6p - fll.7p are quite simple and their 
explicit dependance on the physical parameters is transparent. Moreover, a careful analysis 
of the constants involved in (11. 6p yields the following more explicit two-sided estimate, 

n— 1 / ^ \ n—1 

|r| <dimH^w<dim^^VF<4(l + ^^^|r|^/("-^M , (1.8) 

in all space dimensions n > 3. It is worth pointing out that the bounds in (II. 8p are propor- 
tional to the "surface area" |r| of F, and not the "volume" \Q\ of Q. This is remarkable; most 
nonlinear equations arising in mathematical physics, involving the Laplacian on bounded 
domains, have the dimension of the attractor of the order of \^\" , for some a > and for 
sufficiently large domains. This property may have profound implications in the prediction 
of weather and climate. The reader is referred to Section 4 where this interesting physical 
observation is further discussed for the balance equations governing the large-scale oceanic 
motion. 

Our paper is organized as follows. In Section [21 we obtain upper bounds (cf. Theo- 
rem 12. 7p for the fractal dimension of the global attractor for equation (II. ip with dynamic 
boundary conditions of the form (II. 4p . In Section [3], we employ the same technique of [3] 
to derive a lower bound for the dimension of the unstable manifold of a constant stationary 
solution u* of (II. ip . (II. 4p . As a consequence, we find a lower bound on the dimension of 
Aw (see Theorem 13. ip . Finally, in the Appendix, we recall some useful results on the so- 
called Wentzell Laplacian, and prove an auxiliary inequality, namely, we derive some kind 
of Sobolev-Lieb-Thirring inequality that is required to prove the upper bound in (II. 6p . 
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2 Upper bounds on the dimension 

We use the standard notation and facts from the dynamic theory of parabohc equations (see, 
for instance, [11], [17], [22]). We denote by and IHIpp' norms on (Q) and 
LP (r) , respectively. In the case p = 2, (-, ■)2 stands for the usual scalar product. The norms 
on (Q) and (T) are indicated by |H|j:^r(Q) and |H|j:^r(r), respectively, for any r > 0. 
The natural phase-space for problem (11. ip . (11.41) is 

:= L^iQ) © L^(r) = {F = Q : / G L^iQ), g G L^T)}, 

for all p G [1, oo], endowed with the norm 

\\F\\l, = ljfix)fdx + ljgix)f^, b>0, (2.1) 

if p G [1, oo), and 

||-^||x°° := ||/||L°°(f7) + ||fi'|U°°(r)- 

In the definition of X^, dx denotes the Lebesgue measure on Q, and dS denotes the natural 
surface measure on F. Moreover, we have |llj . 

XP = LP {Q,dfi), pe [l,oo], 

where the measure dfj = dx\Q © ^|p, on f2, is defined for any measurable set 5 C by 
/i(5) = \Bnn\ + \BnT\. The Dirichlet trace map Tr^ : (H) (T) , defined by 

Tr£) (u) = U|r extends to a linear continuous operator Tvd : (il) — )■ H'^^^^'^ (F) , for all 
r > 1/2, which is onto for 1/2 < r < 3/2. This map also possesses a bounded right inverse 
Tr^^ : /f^-V2 (r) ^ jjr (^Q^ g^p^ ^^^^ rj.^^ {Tr~^ij) = tp, for any ^ G //^"Vs (r). Identifying 

each function v e C (U) with the vector V = G C (H) x C (F), it follows that C (IT) 
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is a dense subspace of X^, for every p G [1, oo), and a closed subspace of X°°. Finally, we can 
also introduce the subspaces of H'^ x if'"^V2 (jp^ ^ 

Yr := I Q e {Q) X H^-'/' (r) : Tm (u) = V^j , 

for every r > 1/2, and note that we have the following dense and compact embeddings 
Yr-^ C Yr2, for any ri > r2 > 1/2. The linear subspace V^. is densely and compactly 
embedded into X^, for any r > 1/2. We emphasize that V^. is not a product space and that, 
due to the boundedness of the trace operator Tr^, Yr is topologically isomorphic to H^' {Q) 
in the obvious way. 

We begin by stating all the hypotheses on / and g that we need. We assume that 
g E L'^ (^2) and the following conditions for f E (]R,]R) hold: 

/' (y) > -Cf, for all y eR, (2.2) 

Vi\yf-Cf<f{y)y<V2\yr + Cf, (2.3) 

for some r^^, > 0, C/ > and p > 2. 

We have the following rigorous notion of weak solution to f ll.ip . fll.4p . with initial condi- 
tion u (0) = uq, as in [22] . 

Definition 2.1 The pair U (t) = (^(^j) is said to be a weak solution if ip (t) = Tr^ (n) for 
almost all t E (0, T) , for any T > 0, and U fulfills 

u eC{[q,t]- x^) n (0, T- Vi) n LP (^] X (o, t)) , (2.4) 

u E HIM 00; m, ^ E HIM 00; L' (r)), 
dtUEL'' (0,T;VD, 
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such that the identity 

{dtU, S)x2 + ly (Vm, Va)2 + (/ (u) - An, a)^ = {g, a)^ , 
holds almost everywhere in (0,T), for all S = (^) G Vi. Moreover, we have, in the space 

U (0) = (^J^ =: Uo, (2.5) 

where u (0) = Uq almost everywhere in Vt, and v (0) = vq almost everywhere in T. Note that 
in this setting, vq need not be the trace of Uq at the boundary. Thus, in this context equation 
is interpreted as an additional parabolic equation, acting now on the boundary V. 

The following result is a direct consequence of results contained in [2^ Section 2]. The 
proof is based on the application of a Galerkin approximation scheme which is not standard 
due to the nature of the boundary conditions (see, also, [5]). 

Theorem 2.2 Let the assumptions of Ii2.2\) . Ii2.3\) be satisfied. For any given initial data 
f/o G the problem ( fi. il) . 1^2. 5\} has a unique weak solution which depends continu- 

ously on the initial data in a Lipschitz way. The following estimate holds: 

t+i 

\\Uml + J {\\Uis)\\l^ + \\uis)\\%^^^)ds (2.6) 
t 

<c{\\Uml^)e-'' + c{l + \\g\\l,^^^), 

for allt>0, where the positive constants c, p are independent of time and initial data. 

As a consequence, problem (11.11) . (II. 4p . (12.51) defines a (nonlinear) continuous semigroup 
St acting on the phase-space X^, 

5t : X^ ^ X^ t > 0, 
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given by 



StUo = U{t). 



Theorem 2.3 Let f satisfy assumptions l[2^) . ^KE), let g e L°° {^l) and T e . Then, 
{St} possesses the connected global attractor Aw, which is a bounded subset of Y2 nX°°. As 
a consequence, the global attractor contains only strong solutions. 

Proof. The existence of an absorbing set in Vi fl (Q) and, hence, the existence of the 
global attractor Aw C Vi follows from [22l Theorem 2.8 and Corollary 3.11]. We will now 
show that the attractor is bounded in X°°, and also in V2. All the calculations below are 
formal. However, they can be rigorously justified by means of the approximation procedures 
devised in [22] and [16] (cf. [5] also). From now on, c will denote a positive constant that is 
independent of time and initial data, which only depends on the other structural parameters 
of the problem, that is, \Q\ , \T\ , r]^, u, b, \\g\\^ and n. Such a constant may vary even from 
line to line. 

Step 1. We will first establish the existence of a bounded absorbing set in X°°. First note 
that by (12. 6p . there is a constant Cq > 0, independent of time and initial data, such that for 
any bounded subset B of X^, 3 r = r (||_B||jj2) > with 



sup||f/(t)||x2 <Co. 



(2.7) 



t>r 



We shall now perform an Alikakos-Moser iteration argument. We multiply (II. ip by \u 
Tk '■= 2'', k > 1, and integrate over Q. We obtain 



r-fc-2 



U, 




ufrl + (/ (u) , l^r""^ «>2 + / ^^-^ u) dx (2.8) 



u / dnu\ip\'^'' ipdS + (\u + g,\u\'^'' "^u)^. 



Jr 
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Similarly, we multiply (11. 4p by I'j/'p 'ip/h and integrate over F. We have 



(2.9) 



Adding the equalities (12. Sp . (12. 9p . we deduce 



-^A\\mr,) + {f{u)M'-'n)^ + v I Vu-V{\ur-'u) dx 
M Jn 



{Xu + g^lul''" u)^. 



(2.10) 



A simple manipulation of the third integral in (12.101) . and employing assumption (12. 3p on /, 
we readily get the estimate: 



di 



(ii^iix'j+wikii::S:2+^(2'-i)2 



2-fe 



(2.11) 



<rk{\u + g + Cf,\u\'''' u)^. 



Next, using the fact that |yp ^ < |?/|'"'' + 1, for all /c > 1 and j/ G M, we estimate the last 
term on the right-hand side of (12.111) . 



{\u + g + CjAur-'u)^<c{\\utl + l), 



(2.12) 



for some positive constant c that depends on A and the L°°-norm of but is independent of 
k. On the other hand, it follows from Gagliardo-Nirenberg inequality, and Young's inequality 
for e e (0, 1) , that 



II ' 



(2.13) 



12 



which imphes 



|2 ^ 1 ~ ^ II ||2 -n/2-1 II ||2 



Note that the estimate fl2.13p remains vahd if one replaces the (fi) and (i7)-norms by 
the (r) and (r)-norms, respectively, and n by n — 1, respectively. Setting v = ImI*^*"^ in 
the above inequality, noting that (2'^ — l) 2^"^ > 2, for each k, and the fact that Ttd maps 
(Q) boundedly into (F), we can estimate the gradient term in (12.111) in terms of 



e 



ii,r 



(see, e.g., [351 Chapter 5]). This estimate together with (12. lip . (I2.12p yield 



mr^r,) + c (u^—^ - n) (ikii:: + ii^ii::,r) (2.14) 



for all > 1, where c > is independent of k. 

We shall now make use of an iterative argument to deduce the existence of a bounded 
absorbing set in X*"*, for all k > 1. Thus, noting that rk < (r^)"'''^'''^, then choosing e = 6/rk 
with small 6 = 6 (u) > such that 



1 

r-k] > rk, 



and setting 

yk it) := ^ \u (t, ■)r dx+jjij (t, ■)r ^ = mr^r, , (2.15) 

from (I2.14P we derive the following estimate: 

dtyk (t) + cnyk (t) < c (r,)"/'+^ CVfc.i (t) + 1)' . (2.16) 



13 



Let us now take two positive constants t, /i such that t — ^/r^ > 0, for all k > 1. Their 
precise values will be chosen later. We claim that 

yk{t)<Mk{t,fi):=c{nr^''-'{ sup yk^i{s) + lf (2.17) 

s>t-fi/rk 

holds for y^, defined by fl2.15p and k > 1. To this end, let ( (s) be a positive function 
C : M+ ^ [0,1] such that ( (s) = for s G [0,t-/i/rfc], C (s) = 1 if s G [t,+oo) and 
\d(/ds\ < Cvk, if s G (t — fi/rk,t). We define (s) = ( (s) y^ (s) and notice that 

d d 

— Zk (s) < CTkZk (s) + C (s) — J^A: (s) . 

ds ds 

Combining this estimate with (12.161) and noticing that Zk < 3^fc, we deduce the following 
estimate for Z^'- 

^Zk (s) + cr-kZk (s) < Mk (t, /i) , for all s E [t - /i/r^, +oo) . (2.18) 

Integrating (I2.18P with respect to s from t — fi/rk to t and taking into account the fact that 
Zk {t — fJ'/rk) = 0, we obtain that yk (t) = Zk (t) < Mk {t,fi) (l — e"*"^), which proves the 
claim f l2Tfl) . 

Let now r' > r > be given with r as in (12. 7p . and define /i = 2(r' — r), to = and 
tfc = tk-i — fJ^/rk, k > 1. Using (I2.17p . we have 

sup yk (t) < c (rfc)"/'+^ (sup yk^i (s) + 1)2, A; > 1. (2.19) 
*>*fc-i «>tfe 

Note that from (12. 7p . we have (sup^>j^^^ (s) + l) < Co + 1 =: C. Thus, we can iterate in 



14 



fl2.19p with respect to A; > 1 and obtain that 



sup yk (t) < c{rk) 
t>tk-i 



< ^Ak2Bkn/2+l ^^yk 



cinf'^' (C) 



where 



and 



Therefore, 



Afc := 1 + 2 + 2^ + ... + 2^= < 2^= - 



i=l 



Bk := k + 2 {k - 1) + 2^k - 2) + ... + 2'' < 2''J2 



i=l 



I 

2}' 



SUp3^,.(t) < sup yk{t) < c^k2Bk(n/2+l)-^k ^ 
t>to i>ifc-i 



(2.20) 



(2.21) 



(2.22) 



Since the series in fl2.20p and f l2.2ip are convergent, we can take the r^-root on both sides of 
(1^^ and let A; +CX). We deduce 



sup ||t/(t)||xoo < hm sup(3^fe(t))'/''^ < Ci, 
t>to=T' t>to 



(2.23) 



for some positive constant Ci independent of t, k, U and initial data. 

Step 2. We claim that there is a positive constant C2, independent of time and initial data, 
and there exists r" > such that 



|f/(t)|lv, <C2, for all t>T . 



(2.24) 
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Before we prove fl2.24p . let us recall the following estimate (see [221 Theorems 3.5, 3.10]): 



snp(\\Uit)\\l^ + \\dtu {t)\\l + ]-\\dti^ 

t>TO \ 



2 

2,r 



(2.25) 



t+i 



+ sup 

t>TO Jt 



l'9tM(s)||j:^i(f^) ds 



for some positive constant C3 that is independent of time and the initial data. In order to 
deduce fl^:^ from (1^:^ and (OS]), we need to differentiate dH]) and ([H]) with respect to 
time. This yields 

dfu = lyAdtu - /' (u) dtu + Xdtu, {d^tp + ubdn (dtu)) = 0. (2.26) 

Then, we multiply the first equation of (12.26P by d^u{t) and integrate over VL, using the 
boundary condition of (12.261) . After standard transformations, we obtain 

[W^dtu ml) + pin {t)\\l + \ pli; {t)\\l^ 
= -((/' {u{t))-X)dtu{t),d',u{t))^. 

Using Holder and Young inequalities, we have 



|(||V9,«(t)||^) + ||9^(t)||^ + |||9,2V^(t)||^_^ 
<c( f {u{t))dtu{t) \\\dtu 



<Q{\\u 



00) \\dtU{L)\\^, 



for some positive nondecreasing function Q that depends only on / and c. This estimate 
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yields, owing to (E^l, fl2:25|) . 

j^\\vd,uml<c. 

Then, we can apply the so-called uniform Gronwall's lemma (see, e.g., |l6l Chapter III, 
Lemma 1.1]) to find a time Ti > 1, depending on Tq and r, such that 

\\VdtU {t)\\l < c, for all t > n. (2.27) 

Therefore, (12.271) and (12.251) allow us to deduce from (II. ip and (11.41) . via standard elliptic 
regularity, the following estimate 

||w(t)||^2(^)<c, Vt>ri. (2.28) 

Summing up, we conclude by observing that (I2.24p follows from (I2.28P and the boundedness 
of the trace map Tr£) : [Vl) — t- H^^^ (F). This completes the proof of the theorem. ■ 

Remark 2.4 The proof of Theorem \2.3\ shows how to get an absorbing set in V2. Be- 
cause of this, we can also prove the existence of a global attractor for the dynamical system 
({5,h>o,Vi). 

Theorem 2.5 If Q is a bounded C°° -domain, and f,g are functions, then the global 
attractor Aw is a bounded subset ofYk, for every k > 1. In particular, if U E Aw then 
ueC°° (H) . 

The proof of this result is standard and follows by successive time differentiation of the 
equations in (I2.26P and an induction argument. We omit the details. 

To prove the finite dimensionality of the global attractor Aw-, we can proceed in two 
different ways. One way is to establish the existence of a more refined object called expo- 
nential attractor £wi whose existence proof is often based on proper forms of the so-called 
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squeezing/ smoothing property for the differences of solutions. This can be done by assuming 
smoother nonhnearities, i.e., / G (M) (see, e.g., [ISllTT]). This has been carried out in [16j . 
and references therein, for a system of reaction-diffusion equations with dynamic boundary 
conditions of the form (11.41) . without relating the attractor dimension to the physical pa- 
rameters of the problem. However, since we wish to find explicit estimates of fractal or/and 
Hausdorff dimension of Aw, we shall employ the classical machinery for proving the finite 
dimensionality of the global attractor Aw- This is based on the so-called volume contrac- 
tion arguments and requires the associated solution semigroup St to be (uniformly quasi-) 
differentiable with respect to the initial data, at least on the attractor (see, e.g., [3]). 

We give without proof the following result, which follows as a consequence of the bound- 
edness of Aw into V2 fl X~. 

Proposition 2.6 Provided that f E (M) satisfies the conditions li2.S\) and Ii2. 3\) . the flow 
St generated by the reaction- diffusion equation (iXjy and dynamic boundary condition il.4\) 
is uniformly differentiable on Aw-, with differential 

L (t, U (t)) : = {^^^ = (^^ e X\ (2.29) 

where V is the unique solution to 

dtv = uAv - f' {u{t))v + Xv, ((9i¥? + z/69nf)|r = 0, (2.30) 

v{o) = e. 

Furthermore, L {t, U {t)) is compact for all t > 0. 

The main result of this section is 
Theorem 2.7 Let the assumptions of Proposition be satisfied. The fractal dimension of 
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Aw admits the estimate 



dimp Aw < Co 1 + 



Cf + X 



Cw i^, r) u 



n-1 



for n > 2 



{2.31] 



and 



dimp Aw < Co ( 1 + 



Cf + X 



1/2 



for n = 1, 



(2.32) 



where Cq depends on the shape of Q only. The positive constants Cw^Cd depend only on n, 
Q, r, b and are given in the Appendix. 



Proof. In order to deduce fl2.3ip - fl2.32p . it is sufficient (see, e.g., ^ Chapter III, Definition 
4.1]) to estimate the j-trace of the operator 



L{t,U{t)) 



^ uA- f {u{t)) + XI ^ 



V 



-budr, 







We have 



Trace (L (t,f/(t))Q, 



'X2 



x:(L(t,f/(t))^^,^^.>^ 
j=i 

m m 

i=l i=l 
m m 



1=1 



i=l 



where the set of vector-valued functions • G fl Vi is an orthonormal basis in Q^^. 



Since the family (p^ is orthonormal in Qm^^, using assumption (12.20 on / (i.e., /' [y] > — c/, 
for all ?/ G M), we find 



Trace {I. {t,U) < J] || V</^, H' + (c/ + A) 



m. 



1=1 
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Let n>2. From fl5.12p (see Appendix, Proposition 15. 5p . we obtain 

Trace (L (t, U) Q„) < -vciCw T) m^+^ + {ciuCw T) + c/ + A) m 
=: p (m) . 

The function p (y) is concave. The root of the equation p{d) = is 



Thus, we can apply [H Corollary 4.2 and Remark 4.1] to deduce that dim p Aw < d*, from 
which fl2.3ip follows. The case n = 1 is similar. ■ 

Remark 2.8 Concerning the reaction- diffusion equation U.l\) . we can also handle dynamic 
boundary conditions that involve surface diffusion: 



where a > and Ar is the Laplace- Beltrami operator on T. Our method of establishing upper 
bounds, comparable to the bounds li2.31\) - f2.3S\) . for the dimension of the global attractor can 
be also extended to this case as well. The details will appear elsewhere. 

3 Lower bounds on the dimension 

Lower bounds on the dimension of the global attractor are usually based on the observation 
that the unstable manifold of any equilibrium of the system is always contained in the global 
attractor (see, e.g., [3j). Thus, a lower bound on the dimension of the attractor Aw can 
be found by analyzing the dimension of an unstable manifold associated with a constant 
equilibrium Z for (11.11) . (11.41) . We begin by assuming that g is constant, for the sake of 




dtu — aAru + bi'da4> = 0, on F, 



(2.33) 
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simplicity. Steady-state solutions of (11. ip . (11. 4p satisfy 



Lq {u) := uAu — f (u) + Xu — g = 0, {dnu)\j- = 0. 

We seek a solution of this system U = {xr^^u)) ^ which coincides with a constant 
vector Z = c = (^), c is a constant. Such a stationary solution satisfies the equation 
Lq (z) := —f (z) + Xz — g = 0. Since 

f{y)y> Vi \yf - Cf, for p > 2, 

we have Lq (z) z < —rji \zf + Cf, for some positive constants rji,Cf. Thus, Lq (z) z < 
on the interval Ir = {—R, R) , if R is large enough. It follows that Lq (z) = has at least 
one solution Z = Z {X) (see, e.g., |6, Chapter III]). By the implicit function theorem, this 
solution is of order 1/A for sufficiently large A. 

Now fix this solution. In order to find a lower bound on the dimension of the global 
attractor Aw, it suffices to establish a lower bound for dimE^ (Z) , where (Z) is an 
invariant subspace of L (Z) , which corresponds to 



L{Z)W 



uAw — f {z) w + Xw 
—hvdr.w 



with a (L (Z)) C {C : C > 0}. We note that (L (Z) , D (L (Z))) is self-adjoint on with 
spectrum contained in (— cxd, c/ + A] . 

The main result of this section is the following. 

Theorem 3.1 Let f & C"^ (M) satisfy assumptions Ii2.2\) - I[27^) . There exist a positive con- 
stant cq, depending on f, g and the shape ofQ, independent of X, u, h, , |r| , such that 



dimp Aw > dimn Aw > dimi?+ {Z) > cq 
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A ^ 



Cw i^, r) u 



forn > 2. In one space dimension, the same estimate is valid with Cy/ replaced by Co and 
n — 1, replaced by 1/2, respectively. 

Proof. For a fixed constant solution Z = c oi Lq{z) =0 and sufficiently large A > 1, we 
have X {^) ■= -/' (z) + \ > 0. 

Let be an orthonormal basis in consisting of eigenfunctions of the 

Wentzell Laplacian Aw (see Appendix, Theorem 15. 3p . 

Aw^, = Aj^,, j e No, (p, G D (Aw) n C (H) (3.1) 

such that 

= Ao < Ai < A2 < - < Aj < Aj+i < .... 

We shall seek for eigenvectors Wj = {rp^^^^.^.-j) G X^, of the form Wj (x) = (fj ix)pj, pj E M, 
satisfying equation 



L (Z) W, = (M„ W, eD{L (Z)) := D (A^) • 



(3.2) 



Note that for Wj E D (L (Z)) C Vi, the trace of Wj makes sense as an element of H^^"^ (T). 
Substituting such Wj into (13.21) . taking into account (13. ip and the fact that 

L (Z) W, = -uAwWj + UxW„ U^W, := "^'^ , 

we obtain the equation 



[-iyAjI + Ux)pj = CjPj, n, 



^X(A) 0^ 



v 







(3.3) 
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A nonzero pj exists if C = Cj is a root of the equation 



det {-uAjl + Ha - a) = 0, C > 0. (3.4) 

When z/ = 0, this equation has at least one root C > since x = xW > ^ (iii fact, 
C = xW)- Therefore, there exists 6 > such that when uAj < 6, the equation (13 ■4p has 
a root (j = (j (z/) with (j > 0. Therefore, to any such root (j, we can assign a nontrivial 
Pj, which is a solution of (I3.3p . and thus an eigenvector Wj = {rp^J^), Wj = 'fjPj- Let us 
now compute how many j's satisfy the inequality uAj < 6. The asymptotic behavior of Aj 
is Aj ~ Cpi/ r) j^/'^"'"^^ as j — > oo (see. Appendix, Theorem 15 .4^ . The last inequality 
certainly holds when 

1 < J < ci5"-i (Ch/z/)'"" = C2 (y^) , forn > 2 

and 

/ 1 \ 1/2 

1 < J < c,5'^' {Cj,v)-^l^ = C2 i-^j , forn = 1. 

The positive constants Ci, C2 depend on A. In order to get more explicit estimates for ci, C2, 
it is left to remark that equation (13 ■4p may be rewritten in the form 

det {-uAjX-^I + X-^Ux - Cj) = 

with d = X^^Cj and to observe that a solution of this equation clearly exists if uAjX^^ < S, 
for sufficiently large A and small 6. Employing the asymptotic formula for A^ once again, we 
find 




, for n >2 
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and 




A 



) 



1/2 



1/2 



for n = 1. 



Coy 



It follows that 



dimE+(Z(A)) > 4 



( 



A 



) 



for n> 2 



and 



dimi?. 



(Z(A))>c;ai/2 {Cdv) 



1/2 



in one space dimension. The proof is complete. ■ 

4 Concluding remarks 

In the textbook literature on theoretical geophysics, it was traditional to use a Robin bound- 
ary condition with a nonlinear heat equation to describe temperature variations at the upper 
surface of the ocean [22]. But it was recognized that this was not always the physically 
correct boundary condition [IQ]. Among its applicability to a wide range of phenomena, 
including phase-transitions in fluids, and so on [161 H2] , the reaction-diffusion system (11. ip - 
(11. 4p has important applications in climatology and is essentially used to determine large 
and rapid temperature changes in the ocean's surface as a response to changes into deep 
water formations [ID]. In this paper, we provide explicit bounds for the dimension of the at- 
tractor for this system and study the effect of the dynamic term b^^dfU, representing change 
in thermal energy in an infinitesimal layer near the surface. Unlike the previous results, 
the dimension of the attr actor is proportional to the surface area |r| , for large domains Q 
and fixed parameters u, A and b. Moreover, all the constants involved in our estimates are 
given in an explicit form. We also observe that in the case without h~^dtu in (II. 4p . i.e., 
h = +00, the dimension of the attractor is much larger (and proportional to the volume 
of Q) than the dimension of the global attractor for the same system when < b ^ +oo. 
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Thus, we observe that the addition of the dynamic term b~^dtu, b > drastically changes 
the situation. This is a remarkable fact that can have a profound effect onto the long-term 
dynamics of other systems that are subject to dynamic boundary conditions of this form. 
We will investigate these effects for other systems, such as the Benard problem for nonlinear 
heat conduction, in forthcoming papers. Finally, we note that it is also possible to extend the 
results of this paper to the case when the boundary F consists of two disjoint open subsets 
Fi and F2, each FjVFj is a S'-null subset of F and F = Fi U F2 with Fi ^ F, such that u 
satisfies a Dirichlet boundary condition on Fi and a dynamic boundary condition on F2. We 
will come back to this issue in a forthcoming article. 

5 Appendix 

In this section, we shall recall several important results concerning a certain realization of 
L = uA with the Wentzell boundary condition (11. 3p . We have the following. 

Theorem 5.1 Let Q be a bounded open set o/M" with Lipschitz boundary F. Assume that 
b > and < q E L°° (Q). Define the operator A^ on X^, by 



Then, (A^y, -D (A^y)) is self-adjoint on X^. Moreover, the resolvent operator (/ + A^y) G 
C (X^) is compact. 

We refer the reader to [U UHl HH] for an extensive survey of recent results concerning the 
"Wentzell" Laplacian Aw- 




(5.1) 



with 




(5.2) 
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The eigenvalue problem associated with the operator A\y is given by Awy^ = ^<f', this 
leads to the following spectral problem for the perturbed Laplacian 



— Aip + q (x) If = Aip in Q, (5.3) 

with a boundary condition that depends on the eigenvalue A explicitly: 

bd^ip = Kip on r. (5.4) 

Such a function ip will be called an eigenfunction associated with A and the set of all eigen- 
values A of fl5.3p - fl5.4l) will be denoted by A^y. Let and A^y.,-, j E J, denote all the 
eigenf unctions and eigenvalues of f l5.3p -( l5^ . We have the following (see, e.g., [21 Ii5]). 

Theorem 5.2 Let q > with J q {x) dx > 0. Then, there exists a sequence of numbers 

n 

< Aw,i < Aw,2 < - < Kw,j < ^wj+i < -, (5.5) 

converging to +oo, with the following properties: 

(a) The spectrum of Aw is given by 

a (Aw) = {Aw,j}jm^ 

and each number A\yj, j G N, is an eigenvalue for A^ of finite multiplicity. 

(b) There exists a countable family of orthonormal eigenf unctions for Aw which spans ^ . 



26 



More precisely, there exists a collection of functions with the following properties: 

ipj e D {Aw) and Aw(Pj = ^wj^j, j e (5.6) 
{v3^Vk)^2 = ^jk, j,k e 

= (Blin.span {^^j] (orthogonal direct sum). 

(c) IfV is Lipschitz, then every eigenfunction (fj e Vi, and in fact (fj e C(Q) fl C°°(Q); 
for every j. IfV is of class C^, then every eigenfunction (pj e Vi fl (fi) , for every j. 

(d) The following min-max principle holds: 

^w,3 = min max Rw {(fi, (fi) , j e N, 

dim Yj =j 

where the Rayleigh quotient Rw, for the perturbed Wentzell operators, is given by 

Concerning the case ? = 0, we have the following. 
Theorem 5.3 Let q = 0. Then, there exists a sequence of numbers 

= A^y^^o < ^w,i < ^W,2 < ■■■ < ^wj < ^W,j+l < 

converging to +oo, with the following properties: 
(a) The spectrum of Aw is given by 

O- {Aw) = {^wA j&Vj{Q} ' 

and each number A^kj, j G Nq = NU {0} , is an eigenvalue for Aw of finite multiplicity. The 
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(5.7) 



(5.8) 



eigenvalue Avi/,o is simple and its associated eigenfunction is of constant sign. 

(b) There exists a countable family of orthonormal eig enf unctions for Ayy which spans X^. 
More precisely, the same conclusion (b) of Theorem \5.S\ holds in this case as well. Finally, 
both conclusions (c) and (d) in Theorem \5.2\ hold in the case g = as well. 

The asymptotic behavior of the eigenvalues Awj, as j — t- oo, was estabhshed in [131 
[H]. We refer the reader to \l8l for more details about the Wentzell Laplacian and other 
generalizations. Let J = No or N, according to whether g = or g > respectively. Set 

Cd (fi) := -^^^ and Cs (T) 



Here f„ denotes the volume of the unit ball in M", and we recall that \Q\ stands for the 
ra-dimensional Euclidean volume of Q, while |r| stands for the usual {n — l)-dimensional 
Lebesgue surface measure on T. 

We summarize these results in the following. 

Theorem 5.4 The eigenvalue sequence {Ai^/,,}^.^^ of the (un)perturbed Wentzell Laplacian 
/S.W satisfies: 

(i) For n > 2, we have 

AH.,, = C^(fi,r)jV(-i)+o(jV("-i)), as J ^+oo, (5.9) 

for some 



6C5(r) r2-i/("-i),ll , forn>3, 
Cw{n,T)e{ , ^ ^ ' " (5.10) 



CD{n)Cs{r) 
2{6-iCo{n)+Cs{r))' 



mm{Cnin),bCsiT)} 



, for n = 2. 
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(a) For n = 1, we have 



^w,j = Cd {^) f + o (f) , asj^ +00. (5.11) 

The following version of the Lieb-Thirring inequality is essential. 

Proposition 5.5 Let ujj, 1 < j < m, be a finite family ofWi, which is orthonormal in X^. 
We have 

m 

^ WVujWl > ciCw{^,T) (^m^+^ -m) . (5.12) 

2=1 

The constant ci > depends only on n and the shape of VL, and is independent of the size of 
fi, r, of m, and of the Uj 's. 

Proof. Let Bw := Aw + Cw T) / and let D {Bw) = D (A^/). By Theorems EH ESI 
By/ is a linear positive unbounded self-adjoint operator on X^, such that B^ is compact. 
Thus, we can apply the abstract result of j46, Chapter VI, Lemma 2.1] to deduce that 

m m 

^{\\Vuj\\l + Cw\\uJj\\^^2) = ^{BwUJj,ujj)^2 (5.13) 

1 = 1 4 = 1 

> Avi/,1 {Bw) + Avi/,2 {Bw) + ... + J^w,m {Bw) 

> Cw (1'/^"-') + 2i/("-i) + ... + mi/("~i)) 

> CoCwm~~^'^ , 

since, by (l5.9p -f l5rTT]) . Awj {Bw) > Cw (^5 T) for all j, and some positive constant 

Co (indeed, we have Awj {Bw) = Awj {Aw) + Cw)- Thus, the proof of f l5.12p follows imme- 
diately from fl5.13p . ■ 
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